Spontaneous Symmetry Breaking at the Fluctuating Level 
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Phase transitions not allowed in equilibrium steady states may happen however at the fluctuat- 
ing level. We observe for the first time this striking and general phenomenon measuring current 
fluctuations in an isolated diffusive system. While small fluctuations result from the sum of weakly- 
correlated local events, for currents above a critical threshold the system self-organizes into a coher- 
ent traveling wave which facilitates the current deviation by gathering energy in a localized packet, 
thus breaking translation invariance. This results in Gaussian statistics for small fluctuations but 
non-Gaussian tails above the critical current. Our observations, which agree with predictions derived 
from hydrodynamic fluctuation theory, strongly suggest that rare events are generically associated 
with coherent, self-organized patterns which enhance their probability. 
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Fluctuations arise in most physical phenomena, and 
their study has proven once and again to be a fruitful 
endeavour. The first example is probably Einstein's de- 
termination of molecular scales on the basis of the fluc- 
tuating behavior of a mesoscopic particle inmersed in a 
fluid [1], which opened the door to an experimental ver- 
ification of the molecular hypothesis. Other examples 
range from the role of fluctuations to understand critical 
phenomena beyond mean field theories, to the study of 
fluctuations of spacetime correlations in glasses and other 
disordered materials, which has revealed the universal ex- 
istence of dynamic heterogeneities in these systems [2 . 
In all cases the statistics of fluctuations encodes essential 
information to understand the physics of the system of 
interest. Even further, fluctuations reflect the symme- 
tries of the microscopic world at the macroscale. This 
is the case for instance of the Gallavotti-Cohen fluctua- 
tion theorem [SHS] or the recently introduced isometric 
fluctuation relation |6], which express the subtle but en- 
during consequences of microscopic time-reversibility at 
the macroscopic level. Special attention is due to large 
fluctuations which, though rare, play a dominant role as 
they drastically affect the system behavior. 

The study of fluctuating behavior provides an alter- 
native way to derive thermodynamic potentials from 
which to calculate the properties of a system, a path 
complementary to the usual ensemble approach. This 
can be extended to systems far from equilibrium |2UTT], 
where no bottom-up approach exists yet connecting mi- 
croscopic dynamics with macroscopic properties. The 
large-deviation function (LDF) controlling the fluctua- 
tions of the relevant macroscopic observables plays in 
nonequilibrium systems a role akin to the equilibrium 
free energy, and reflects the phenomenology typical of 
nonequilibrium physics (e.g. non-local behavior result- 
ing in long-range correlations [3 Ej). Hydrodynamic 
fluctuation theory (HFT), which studies dynamic fluc- 
tuations in diffusive media [8 11 , offers predictions for 



both the LDF and the optimal path in phase space re- 
sponsible of a given fluctuation, which can be in general 
time-dependent [9l[T0]. However, it has been shown that 
this optimal path is in fact time-independent in a broad 
regime [HI [13] . This scenario eventually breaks down for 
large fluctuations via a dynamic phase transition at the 
fluctuating level involving a symmetry breaking [TOl [H] . 

In this paper we report compelling evidences of this 
phenomenon in a paradigmatic model of transport in 
one dimension (ID), where we study fluctuations of the 
time- aver aged current. We find that small current fluc- 
tuations result indeed from the sum of weakly-correlated 
local random events in the density field, thus giving rise 
to Gaussian statistics as dictated by the central limit the- 
orem, see Fig. [l]a. However, for large enough currents, 
the system self-organizes into a coherent traveling wave 
which facilitates this rare event by accumulating energy 
in a localized packet, see Fig. [l]b, with a critical current 
Qc separating both regimes. This phenomenon, predicted 
by HFT [TOl [14! , is most striking for this model as it hap- 
pens in an isolated equilibrium system in the absence of 
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FIG. 1. (Color online) Typical evolution of the energy field 
for different current fluctuations in the ID KMP model on 
a ring, (a) Small current fluctuations result from weakly- 
correlated local events, (b) However, for \q\ > qc the system 
facilitates this unlikely deviation by forming a traveling wave. 
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FIG. 2. (Color online) Main: Measured /i(A) for the ID KMP 
model on the ring and increasing values of iV, together with 
the HFT prediction and the Gaussian approximation. Inset: 
/i(A) — /ifiat(A) for the same N. Data converge to the HFT 
prediction as N increases. 



any external field, breaking spontaneously a symmetry 
in ID. This is an example of the general observation that 
symmetry-breaking instabilities forbidden in equilibrium 
steady states can however happen at the fluctuating level 
or in nonequilibrium settings [15 . Such instabilities may 
help explaining puzzling asymmetries in nature [15 , from 
the dominance of left-handed chiral molecules in biology 
to the matter-antimatter asymmetry in cosmology. 

Our model system is the paradigmatic ID Kipnis- 
Marchioro-Presutti (KMP) model of transport on a ring 
[16]. This is a general model of transport which repre- 
sents at a coarse-grained scale the physics of many quasi- 
ID systems of theoretical and technological interest char- 
acterized by a single locally-conserved field which diffuses 
across space. In this sense our results are of great gen- 
erality and may have important implications in actual 
experiments. Moreover, this model acts as a benchmark 
to test theoretical advances in nonequilibrium physics 
H [IQI m [H [17]. The model is defined on a ID lattice 
of N sites with periodic boundary conditions. Each site 
z G [1, A^] is characterized by an energy pi > 0, and mod- 
els an oscillator which is mechanically uncoupled from its 
nearest neighbors but interacts stochastically with them 
via a random energy redistribution process which con- 
serves total energy po = Ylf^iPi- We are interested in 
the statistics of the total current q flowing through the 
system, averaged over a long diffusive time r. For r ^ oo 
this time average converges toward the ensemble average 
(g), which is of course zero because the system is iso- 
lated and in equilibrium. However, for long but finite 
r we may still observe fluctuations q / (g), and their 
probability P^- {q) obeys a large deviation principle in this 
limit [18 , Pr(<7) ^ exp[+rA/'G((7)]. This means that the 
probability of observing a current fluctuation decays ex- 
ponentially as both r and N increase, at a rate given by 
the current LDF G{q) < 0, with G{{q)) = 0. For KMP 



model a singularity has been shown to exist in G{q) [10] 
whose details we uncover here. 

In order to study in depth current statistics, we per- 
formed extensive simulations of the ID KMP model with 
Po = 1 using an advanced Monte Carlo method which al- 
lows to explore the tails of the current LDF [191 [20]. This 
method implies a modification of the stochastic dynamics 
so that the rare events responsible of a large current fluc- 
tuation are no longer rare, and requires the simulation 
of multiple clones of the system [19]. In this work we 
used M = 10^ clones. The method yields the Legendre 
transform of the current LDF, /i(A) = m8iXq[G{q) + Ag], 
with A a parameter conjugated to the current, and Fig. 
[2] shows simulation results for /i(A) and increasing values 
of N. As shown below, HFT predicts Gaussian current 
statistics for \q\ < qc = tt -see eq. (|3|, corresponding 
to quadratic behavior in /i(A) = /iflat(A) = A^/2 up to 
a critical |Ac| = tt. This is fully confirmed in Fig. |2] as 
increases, meaning that small and intermediate cur- 
rent fluctuations have their origin in the superposition 
of weakly-correlated local events, giving rise to Gaus- 
sian statistics as dictated by the central limit theorem. 
However, for fluctuations above the critical threshold, 
|A| > Ac, deviations from this simple quadratic form are 
apparent, signalling the onset of a phase transition. In 
fact, as N increases a clear convergence toward the HFT 
prediction is observed, with very good results already 
for N = 32. Strong finite size effects associated with 
the finite population of clones M prevent us from reach- 
ing larger system sizes [21 , but = 32 is already dose 
enough to the asymptotic hydrodynamic behavior. Still, 
small corrections to the HFT predictions are observed 
which quickly decrease with see inset to Fig. [2] 

The phase transition is most evident at the configura- 
tional level, so we measured the average energy profile 
associated to a given current fluctuation [20], see Fig. 
|3| Due to the system periodicity, and in order not to 
blur away the possible structure, we performed profile 
averages around the instantaneous center of mass. For 
that, we consider the system as a ID ring embeded in 
two-dimensional space, and compute the angular posi- 
tion of the center of mass, shifting it to the origin before 
averaging. Notice that this procedure yields a spurious 
weak structure in the subcritical region, equivalent to av- 
eraging random profiles around their (random) center of 
mass. Such spurious profile is of course independent of 

and can be easily substracted. On the other hand, su- 
percritical profiles exhibit a much pronounced structure 
resulting from the appearance of a traveling wave, see 
Fig. [ijb. Top panel in Fig. [3] shows the measured profile 
ujx{x) for different A > Ac and varying N. Again, fast 
convergence toward the HFT result is observed, with ex- 
cellent agreement for A" = 32 in all cases. Bottom panel 
in Fig. |3] shows the measured profiles for A" = 32 and 
different A, which closely resembles the HFT scenario, 
see inset. We also measured the average velocity associ- 
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FIG. 3. (Color online) Top: Supercritical profiles for different 
A and varying and HFT predictions. Bottom: Measured 
profiles as a function of A for N = 32. Inset: HFT prediction 
for the optimal profile ujq (x) . Profiles are flat up to the critical 
current, beyond which a nonlinear wave pattern develops. 



ated to a given current fluctuation by fitting the motion 
of the center of mass during small time intervals At to 
a ballistic law, xcuit + At) — xcuit) = vAt^ see e.g. 
Fig. [T]b, and making statistics for the measured veloc- 
ity. Fig. [4] shows the mean velocity for At = 100 Monte 
Carlo steps as a function of A for increasing values of 
A^, and the agreement with HFT is again very good al- 
ready for = 32 (other values of At yield equally good 
results). Notice that for subcritical current fluctuations 
the velocity is simply proportional to the current, while 
above the critical line the relation becomes nonlinear. 

To understand this behavior, note that the KMP 
model belongs to a large class of diffusive systems which 
evolve in time according to a rescaled continuity equation 



dtp = da:(^D[p]da 



(1) 



Here p{x^t) is the density field, with x G [0, 1], j{x^t) = 
— {D[p]dxP -\- £,) is the fluctuating current, and D[p] is 
the diffusivity (a functional of the density profile in gen- 
eral). The (conserved) noise term ^(x,t), which accounts 
for microscopic random fluctuations at the mesoscopic 
level, is Gaussian and white with {^{x^t)) = and 
{^{x,t)^{x\t')) = N-^a[p]S{x - x')5(t - t'), being 



the mobility functional and N the system size. In par- 
ticular, for the KMP model D[p] = 1/2 and cr[p] = p^ 
[16], and we focus here on periodic boundary conditions, 
p(0, t) = p(l, t) and j(0, t) = j(l, t), so the total mass in 
the system is conserved, po = p{x,t)dx. The proba- 
bility of observing a particular history {p(x, t), j (x, t)}Q 
of duration r for the density and current fields can be 
written as a path integral over all possible noise real- 
izations, {^(x,t)}o, weighted by its Gaussian measure, 
and restricted to those realizations compatible with eq. 
([1]) at every point in space and time. This results in 
P({p, j}o) ^ 6^P{+^^r [p, j]}, with a rate functional de- 
fined by the familiar formula jSHTT] 
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(2) 



with p{x^t) and j{x^t) coupled via the continuity equa- 
tion, dtp -\- dxj = 0. Eq. ([2| expresses the locally 
Gaussian nature of current fluctuations around its av- 
erage behavior, given by Fourier's law. We are inter- 
ested in the fluctuations of the time-averaged current 
q = dt dxj{x^t). The probability of observ- 

ing a given q can be in turn obtained from the path 
integral of P({p, j}o) restricted to histories {p, j}o con- 
sistent with that value of q. This probability scales as 
Pt{q) ^ ex.p[-^rNG{q)]^ and the current LDF G{q) is re- 
lated to Xr[p^j] via a simple saddle-point calculation in 
the long time limit, G{q) = mdiXpj Xr[p^j], such that 
the optimal profiles pq{x^t) and jq{x^t) solution of this 
variational problem are compatible with the constraints 
on po and q and are related via the continuity equation. 
These optimal profiles can be interpreted as the ones the 
system adopts to facilitate a given current fluctuation. 

Small deviations of the empirical current away from its 
ensemble average (q) = typically result from weakly- 
correlated local fluctuations. The average density profile 
associated to these small fluctuations hence corresponds 
still to the flat, stationary one, pq{x, t) = po. In this case, 
the optimal current profile is just jq{x,t) — q and the 
current LDF is simply quadratic, Gflat(^) = —Q^/'^^{po) 
[14^ , resulting in Gaussian current statistics as confirmed 
is simulations, see Fig. [2] A natural question thus con- 
cerns the stability of this flat profile against small per- 
turbations. Bodineau and Derrida have shown [14 that 
the flat profile indeed becomes unstable, in the sense that 
G{q) increases by adding a small time-dependent periodic 
perturbation to the otherwise constant profiles, whenever 
87r'^D{pQ)'^a{po) — q'^a^\po) < 0, where a^^ denotes second 
derivative. This defines a critical current 



27rD{po)V2a{po)/a^^{po) 



(3) 



for the instability to kick in. When this happens, the 
form of the associated relevant perturbation suggests that 
current fluctuations in this regime are sustained by a 
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FIG. 4. (Color online) Velocity measured as a function of A 
for increasing N, and HFT result. 



traveling wave pattern moving at constant velocity v [14] , 
as observed in simulations. We hence write pq{x^t) = 
ujq{x — vt), which results in jq{x, t) = q — vpo-\-vujq{x — vt) 
via the continuity equation. The variational problem for 
G{q) can be now written as 



G{q) = — min 



^ [q - VPO + VOJq{x)] ^ + UJq{x)'^D[uJq]'^ 



fluctuation theory, is most surprising as it happens in 
an equilibrium system in the absence of external flelds, 
breaking spontaneously a symmetry in ID. This illus- 
trates the idea that critical phenomena not allowed in 
equilibrium steady states may however arise in their fluc- 
tuating behavior or under nonequilibrium conditions [15] . 
Remarkably, similar instabilities have been described in 
quantum systems ^ . Our results strongly support that 
the phase transition is continuous as conjectured in [14], 
excluding the possibility of a flrst-order scenario, and sug- 
gest that a traveling wave is in fact the most favorable 
time-dependent proflle in the supercritical regime. This 
observation may greatly simplify general time-dependent 
calculations, but the question remains to whether this is 
the whole story or other, more complex solutions may 
play a dominant role for even larger fluctuations. In any 
case, it seems clear that rare events call in general for 
coherent, self-organized patterns in order to be sustained 
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2a[ujq] 



dx 



(4) 

resulting in the following differential equation for the 
shape of the optimal traveling wave 

[q-vpo^vujq{x)]'^ -ujq{xf D[ujq]'^ = 2(7 [cj^] | Ci +^20;^ (x) | . 

(5) 

This equation yields a 0Jq{x) which is generically a sym- 
metric function with a single minimum Ui = Co'(xi) and 
maximum a;2 = ^{x2) such that |x2— xi| = 1/2. The con- 
stants Ci and C2 can be then related to these extrema, 
which in turn are flxed by the constraints on the total 
mass of the system and the distance between extrema. 



ujD{uj) 



2 Z,{uj) 



duj , (6) 



where Zy{uj) = [{q - vpo + vcu)'^ - 2a{u){Ci + C2Cj)]^/^. 
The optimal wave velocity is given implicitely hy v = 
-Q^["^/4"^ [14, with 



D{uj){uj-po)'' 
a{uj)Zy{uj) 



duj 



In this way, for constant po and q^ we use eqs. ([6|-([7|) 
to compute the proflle extrema and its velocity, and use 
this information to solve eq. ([5| for ujq{x). The resulting 
predictions are fully conflrmed in simulations, see above. 

Our results unambiguously show that an isolated dif- 
fusive system exhibits a phase transition at the fluctua- 
tion level. This phenomenon, captured by hydrodynamic 
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